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ABSTRACT 


A review of the understanding of the combustion mechanisms of solid propellants that the authors have 
built from their work and from the literature is presented. Such an understanding is an important part of 
the process carried out to master the behavior of solid propellants and to obtain desired characteristics 
(with respect to energetic level, burning rate level, sensitivity to pressure and initial temperature, nature 
of emitted combustion products, vulnerability to various aggressions...). 
The propellants and propellant components considered are: 

¢ double-base propellants, based on nitrocellulose and nitroglycerin, 

* active binder, based on an inert polymer (or energetic such as PAG) and a liquid nitrate ester, 

¢ inert binders, such as polybutadiene, 

* ammonium perchlorate, 

e  nitramines, such as HMX, RDX and CL20 (HNIW), 

* composite ammonium perchlorate-inert binder propellants, 

* composite propellants based on a nitramine and an active binder, 

¢ aluminum, with respect to the two previous types of propellants, 


e additives, when appropriate. 


The features of the combustion zone described are: 


¢ In the condensed phase, the thickness of the temperature profile and of the decomposition zone, 
the kinetics of the decomposition, the energy released, the nature of the gases evolved, the surface 
temperature, 


¢ In the gas phase, the type of flame structure (diffusion or kinetically controlled), the possibility of 
staging (such as in double-base propellants), the kinetics of the reaction(s), the energy released, 
the flame temperature (primary and final, when applicable). 


It is concluded that a fairly proper knowledge of the combustion of the various components and 
propellants has been acquired (being now extended to new ingredients, oxidizers or binders). 
Furthermore, based on this knowledge, a first approach modeling description can be achieved. Such a 
description is necessary in accompanying the elaboration of new propellants and in preparing the 
investigation of more complicated regimes such as those of erosive burning and of non-stationary 
response. 


Paper presented at the RTO/VKI Special Course on “Internal Aerodynamics in Solid Rocket Propulsion”, 
held in Rhode-Saint-Genése, Belgium, 27-31 May 2002, and published in RTO-EN-023. 
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INTRODUCTION 


Much work has been devoted in various countries to investigating the combustion mechanisms of solid 
propellants. It is timely to bring together the information obtained by the authors and compared to that of 
the literature on the combustion of the individual components as well as of their combination into 
propellants. This review is about the existing components and propellants: double-base propellants and 
active binders, inert binders, ammonium perchlorate, HMX, RDX and CL20 and the corresponding 
composite propellants, ammonium perchlorate — inert binder (plus possibly aluminum), HMX (or RDX or 
CL20) — active binder. 


The viewpoint adopted here is that of the understanding of the combustion behavior of propellants. 
Therefore as much information as possible is presented about the fundamentals of the processes 
(thermal properties, kinetics in the condensed phase and in the gas phase...), whereas no attempt is made to 
establish a complete catalog of practical results on various propellants with different particle sizes, 
catalysts, variations on the percentage of ingredients. The aim is to give a clear, as conclusive as possible 
picture. This will then be non compatible with a complete discussion of the various, sometimes 
contradictory, mechanisms proposed in the literature. Also precluding such a discussion is the will to 
compare the different components and the corresponding propellants. 


Some space is taken up by physico-chemical modeling. The aim is not so much to give the elements of 
mathematical descriptions which could be used for a priori computations of burning characteristics of 
propellants (to the extent that such computations are possible). The point is more to put to test the 
hypotheses made on the mechanisms of combustion by incorporating them in reasonable models and 
confronting the results thus obtained to experimental data. 


These descriptions can also be viewed, along-side with the data given for each component or propellant, 
as useful for mastering the regimes of combustion which go beyond stationary combustion: that is erosive 


burning and unsteady (under pressure excursions or pressure oscillations) combustion responses. 


Table 1 gives information about the various types of propellants of actual use. 
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Table 1: Performances/Characteristics of Various Propellants 
PROPELLANT | COMPOSITION Pp g/cm? 1,(70/1) theoretical. APPLICATIONS 
(main ingredients) (losses) /CHARACTERISTICS 
Extruded DB Nitrocellulose - Anti-tank rockets and missiles 
Nitroglycerin < 1.66 - AS rockets 
oa - Some tactical missiles (SA) 
Minimum smoke 
Powder cast DB | Nitrocellulose - Anti-tank missiles 
Nitroglycerin < 1.66 - Some tactical missiles (AS) 
Minimum smoke 
AP composite Ex: 88 % AP - 1.72 - Some AS rockets 
12 % HTPB - Some tactical missiles 
Reduced smoke (HCI-H2O) 


AP composite Ex: - AA tactical missiles 
with aluminum | 68 AP - 20 al. - =1.82 - Anti-ship missiles (booster) 
12 HTPB = - Tactical ballistic missiles 
- Strategic ballistic missiles 
- Apogee motors 
- Boosters for space launchers 
(Titan IH, IV, V, Space Shuttle, 
Ariane V, H2A...) 
Smoky (Al,03) 


HMX (RDX) HMX or RDX - < 1.75 - Anti-ship missiles (cruise) 
composite XLDB binder - SA missiles 
Minimum smoke (without AP) 


HMX composite | HMX + AP+ Al - Strategic ballistic missiles 
with aluminum | XLDB binder (upper stages) 
1.8 (Trident, MX...) 
Smoky (A1203) 


Sources: Air et Cosmos n° 1000, May 1984. 
Annales des Mines, Jan-Feb. 1986. 
Aéronautique et Astronautique, n° 138, 1989. 


Double-base propellants (made by the extrusion or powder casting techniques) are used in anti-tank 
rockets or missiles and in some tactical missiles. Their main advantage is that they produce a minimum 
amount of smoke (only from a small amount of additives). 


Composite propellants, based on ammonium perchlorate (AP) without aluminum, generate reduced 
smoke, HCl and H,O vapor will precipitate into droplets in the plume under given temperature and 
humidity conditions. They are used for various tactical missiles. With aluminum, they are widely used 
in missiles and space launchers. They produce alumina smoke, which, in the case of space launchers, 
could be considered in the future to be undesirable (along with HCl). 


Composite propellants based on nitramines and an “active” binder (cross linked polymer with 
nitroglycerin or other liquid nitrate esters) are used more and more. Without aluminum, they are in the 
minimum smoke category and they replace DB propellants. With aluminum, they reach the highest 
specific impulse and density and are used so far for upper stages of strategic missiles. 
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The combustion of the components and then of the various propellants will be seen in the next chapters. 


A few general references about chemical propulsion, solid propellants and combustion can be found at the 
end of the main text, ahead of more specialized references introduced progressively in the following 
chapters. 


COMBUSTION OF DOUBLE-BASE PROPELLANTS AND ACTIVE BINDERS 


1.0 Introduction 


It seems appropriate to consider double-base propellant combustion mechanisms in the first place because 
they correspond to relatively simpler premixed processes which lend themselves to a better understanding 
and because they have been investigated for a long time (starting in the 1950’s) in the US, the then USSR, 
the UK, Japan and France in particular. General background on the combustion of double-base propellants 
can be found in references [1-6]. 


As will be seen, the mechanisms involved in the combustion of double-base propellants will apply as well 
to the active binders ( ~1/3 polymer, = 2/3 nitroglycerin, or other liquid nitrate esters). 


Double-base propellants are made in a number of ways. When they are rolled or extruded, the components 
are nitrocellulose and nitroglycerin, to which some stabilizers such as centralite and plasticizers are added. 
When they are cast, a casting powder (made of nitrocellulose, some nitroglycerin, and the various 
additives) is swelled within the mold by a liquid mixture of nitroglycerin and triacetin. The grain thus 
obtained is then inhibited and used free standing in the motor. The propellant ingredients can also be 
mixed, cast, cross-linked, and the grain case bonded. 


Depending on the relative amounts of nitrocellulose and nitroglycerin, Figure 1, the energetic level of the 
propellant can be increased or, in the usual terminology, its “heat of explosion” or “calorimetric value”, 
that is, the heat evolved in a calorimetric bomb by combustion under an inert atmosphere. One can then 
talk about “cool” and “hot” compositions. 


H, C-ON O, H ONO, 
| | 
C—O O c—cC H 
H l l 
a ST AN 
C ON O, ÇC ON O, Ç 
o | H H H H Lo 
c—cC C—O 
| l x 
H ONO, H, C-ON O, n 
Nitrocellulose 
h 
H, 1 + CH, 
ONO, ONO, ONO, Nitroglycerin 
Example of composition : 
1100 cal/g, 52.25 % NC, 42.75 % NG, 5 % others 


Figure 1: Components of Double Base Propellants. 
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Double-base propellants are used in small and medium sized rockets and thus exposed to varying ambient 
temperatures. The sensitivity of the motor operation to temperature depends upon the propellant burning 
rate sensitivity to both the temperature and the pressure. As can be seen on Fig. 2, the pressure 
exponent, in the usual empirical law v, ~ p”, is around 0.7 and increases to nearly 1 at high pressure. 
Super-rate effects (Fig. 3) are created by the use of additives, most often lead and copper salts combined 
with carbon black. At the end of the super-rate zone, the burning rate falls back to that of the control 
propellant, with the occurrence of a nearly zero pressure exponent zone, a “plateau” effect, or a negative 
exponent zone, a “mesa” effect. These terms are used by analogy with topographical features. A fairly 
complete set of results can be found in reference [6]. It is only in these reduced pressure exponent zones 
that the propellant is used to minimize the motor operation sensitivity to ambient temperature. Due to this 
fact, the study of the combustion of propellants without additives should be conceived only as a first step 
leading to an understanding of modified (that is, with additives) propellants. 


Cool propellant 820 cal/g, 2220 K 
Hot propellant 1100 cal/g, 2835 K 
Very hot propellant 1250 cal/g, 3110 K 
Pure nitoglycerin 1770 cal/g, 3260 K 


> © O 


(Temperatures at equilibrium at 100 atm) 


V, (mm/s) 


100 


50 


10 


0.7 - i 


1100 cal/g 


820 cal/g 


- — jH as 


1 5 10 50 100 500 1000 


Figure 2: Burning Rate vs Pressure. 
(Double Base Propellants) 
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A Cool propellant 
® Hot propellant 
Very hot propellant 
D Cool propellant + 2 parts lead salt 
+ 1 part carbon black 
O) Hot propellant + lead and copper salts 


+ carbon black 
(3') Very hot propellant + lead salt 
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50 100 500 1000 


100 
50 


10 


p (atm) 


Figure 3: Super-Rate Effects. 
(Double Base Propellants) 


2.0 Flame Structure 


From the works mentioned previously it is possible to describe the combustion wave structure of double- 
base propellants, in particular its chemical processes, see Figs. 4 and 5. The various data will be discussed 
and justified later. Gas analysis results are from reference [7]; they refer to mass fractions. 
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T (°C) 


100 atm 
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2000 + 
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Figure 4: Temperature Profiles. 
(Double Base Propellant) 
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SECONDARY (LUMINOUS) 
FLAME 


NO/CO 


© SUPERFICIAL DEGRADATION ZONE 
(FOAM ZONE) 
C - O -//-NOp — 0.467 NO3 


© PRIMARY FLAME (FIZZ ZONE) 


NO% / ALDEHYDES 
NO / CARBON 


INDUCTION (DARK) 
ZONE 


© PREHEATED 
ZONE 


SURFACE 
NO2 0.255 
CH20 0.21 — 
(CHO)p 0.1 -l 
NO 0.14 NO 0.26 
— No 0.04 
0.09 (ee) 0.34 
CO2 0.1 CO2 0.18 
0.06 H20 0.14 
H2 0.007 
HYDROC. 0.04 


CARBON 0.05 


Figures for an 1100 cal/g propellant. Surface and primary flame (at 11 atm) 
mass fractions from gas analysis 


Figure 5: Various Zones of the Combustion of a Double Base Propellant. 


The propellant components pass unaffected through a preheated zone of a few tens of micrometers in a 
few milliseconds and reach a superficial degradation zone (or “foam” zone in the early literature) 
where the temperature becomes high enough for the molecular degradation to take place, initiated by the 
rupture of the C-O-//-NO, bond. Simultaneous recombination occurs so that a mixture of NO,, aldehydes, 
but also NO emerges from the surface and so that the net energy balance of the degradation is exothermic. 
At pressures under about 100 atm, a clearly separated primary flame (“fizz” zone) and a secondary flame 
(“luminous” flame) are observed, the first involving NO -aldehydes reactions and the second probably 
NO-CO reactions. In this pressure range the secondary flame is too far away to have any effect on the 
surface or even to induce a temperature gradient into the primary flame. The burning rate is then entirely 
under the influence of the latter. This corresponds to a burning rate/pressure law with a 0.7 pressure 
exponent (Fig. 1). As the pressure increases, the secondary flame enhances and then merges into the 
primary flame and a transition is observed to a zone with a pressure exponent close to 1. When the 
secondary flame is fully developed, even at pressures for which it does not yet influence the burning rate, 
the final products (N2, CO, CO2, H2O and H2) and the final temperature (2100-3100 K, depending on the 
heat of explosion) are attained. 


Table 2 gives data relative to the various zones of the combustion wave. 
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Table 2: Characteristics of the Combustion Zones (Measured Results from Zenin) 


eooo o 


Preheated zone, um 

(measured/computed) 140/194 50/55 45/35 
Residence time in 100 3 
preheated zone, ms 


Superficial degradation 

zone um 11 3 2 

Residence time in 

superficial zone, ms 0,5 0,2 

Flame thickness, ym 

(measured) 200 75 110 
(secondary flame) 


Double base propellant taged flame (1° flame , dark zone 2 “flare ), under pressures: 
0.7 1.2 1.6 2.1 2.6 MPa. 


HVX uncer 0.15 MPa 


Visualizations with high speed camera, for a double base propellant, showing the secondary flame 
moving closer to the surface, and a comparative picture for HMX which will be evoked further on. 
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3.0 Condensed Phase Processes 


The preheated zone of a regressing propellant is described by the conservation of energy in a coordinate 
(x > 0 in the gas phase) regressing with the surface: 


Pp Vb Cp dT / dx = d( A, dT / dx ) / dx (1) 
in such a way that a temperature profile 
(T-To)/(Ts;- To) = exp( x vs/ dp), dp = Ap/ Pp Cp (2) 


will progress with the surface into the propellant. From measurements up to 100°C and from ignition 
experiments, representative average values are taken as indicated in Table 3. 


Table 3: Values of the Condensed Phase Properties (Double-Base Propellants) 


Pp = 1.6 g/cm’, cp = 0.4 cal/g K, 
Ap = 5.1 10% cal/s cm K 
dp = 0.8 10° cm’/s 


The thickness €cona. of the conduction zone can be taken conventionally as 

T( end of cond. zone ) - Ty = 10° (T;- To) 

cond. = (dy/vp ) In 10° (3) 
As an example, for v, = 10 mm/s 


€cond, = 37pm, a thickness through which the temperature rises from 293 K to about 700 K. The residence 
time through this conduction zone is: 


Teond, = ( dp/ vp) In10° (4) 
about 4 ms in this example; a very short time for a temperature increase of 400 K. 


The superficial degradation zone has its thickness ruled by the conservation of the non degraded 
propellant mass fraction Y,: 


Pp Vo dY, / dx = - pp Ac exp(- E./RT), (5) 
with the decomposition represented by an Arrhenius law. Numerous investigations by thermogravimetry, 


differential scanning calorimetry, on nitrocellulose, nitroglycerin and other nitrate esters, as well as on 
double base propellants, and ignition studies [5] result in: 


Table 4: Values for the Condensed Phase Degradation Kinetics (Double-Base Propellants) 


Decomposition order 0, A, = 1 10” s! E, = 40 kcal/mole 
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The thickness of the degradation layer is related to the fast drop in the degradation rate. When this rate is 
10° that at the surface temperature, the lower limit of the reaction layer is conventionally reached: 


exp[ - & / ( 1 - AT / T;)] = 107 exp( - &), é= E./ RT; 
the temperature drop is then 
AT/T,=1/(1+&/In10 ) (6) 


AT = 100 K for T, = 700 K. Such a temperature drop inserted in Eq. (2) gives an estimation of the reaction 
layer thickness 


Creaction = - ( dp/ Vp) In[ 1 - AT/( T;-T.)] 
or taking into account the magnitude of the reduced activation energy ( & = 30 ). 
Creaction ~ Ccond. Ts / Ge ( Ts a Ty ) (7) 


For the values taken above, at vp = 10 mm/s, erat ~ 2 um, with an associated residence time Treact of 
0.2 ms. 


The summation of Eq. (5) through the degradation layer results in: 
Po Vo( Yps - Ypo ) =- Í Pp Ac exp(- Ec/ R T ) dx (8) 
Vp ~ Ac exp( - Ec/ R Ts) Ereaction = Ac Exp( - Ee) ( dp/ Vb) 1n10°( 1/&)[T;/CTs- To) J 


A more rigorous approach [8] (which is almost identical to a numerical computation [9]) gives, for a zero 
order reaction: 


Vp = (dp / &) Ac exp( - & )/(1-To/T,-Q,/2 c)Ts) (9) 


This equation indicates a relation between surface temperature and burning rate: the mass flow rate pp Vp 
emitted from the surface is the result of the decomposition of the propellant into gases throughout 
the superficial degradation layer, Eq. (8). The higher the burning rate vp, the smaller the residence time 
Treact, ~ 1 / Ve, and the higher is the surface temperature reached to allow for the complete degradation of 
the propellant. 


Traverses with micro-thermocouples (as seen previously the thickness of the combustion wave is of the 
order of tens of um) allow to obtain measurements of the surface temperature [2,3,4,5]. One example is 
given on Fig. 6. The results from various sources are collected (see [5] for references), as burning rate 
versus | / T,, on Fig. 7. Also indicated is the correlation obtained from Eq. (9). Due to the thinness of the 
conduction zone, a few tens of micrometers, fairly large errors and scatter should be expected in the 
measurements of the surface temperature. Nevertheless, some conclusions can be reached (see also 
Refs. [3,10]). The initial degradation of the propellant components is controlled by the breaking of the 
-C-O-//-NO, bond (characterized by the 40 kcal/mole activation energy). This is considered to be a 
temperature sensitive only process, irreversible (therefore not influenced by the pressure level). It is 
noteworthy that the kinetics of the degradation is the same from thermal decomposition (by TG and DSC) 
at about 400 K, to ignition from 400 to 500 K [5] and combustion at temperatures up to 700 K. 
Also important is the conclusion, if one looks at the details of Fig. 7, that the presence of super rate 
producing additives does not affect the condensed phase kinetics. 
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The energetics of the reaction layer is now to be considered. The initial degradation of the propellant, 
taking into account the assumed decomposition of nitroglycerin into 3N02, 2 CH,O and 0.5 (CHO)), 
is thought to give (for the example of the 1100 cal/g propellant) the mass balance 


1 propellant — 0.467 NO, + 0.364 CH,O + 0.1 (CH0); + 0.06 hydrocarbons (10) 
with a corresponding endothermic heat of degradation 
Qu = - 135 cal/g of propellant. 


It is thought [5,9,10] that within the superficial layer the exothermic reaction between NO, and aldehydes 
can start. A plausible mole balance (in order to match various results, in particular the analysis [7] of the 
gases emitted from the surface of regressing propellants) is 


NO, + CH,0 — NO + 1/2 CO + 1/2 CO, + 1/2 H20 + 1/2 H» (11) 
with a corresponding exothermic heat of reaction 
Qno2 = 1040 cal/g of NO2. 


Conservation of the species NO, (in terms of mass fraction Y) is written (no diffusion is taken into 
account) 


Pp Vb dYyno2 /dx= Yno2i Pp Ag exp( = E. /RT ) = Ano? ( p M/RT ) Yno2 exp( a Eno? / RT ) (12) 


if a first order reaction with respect to the molar concentration of NO, is assumed (it will be seen that this 
is probably the case). 


The conservation of energy is written 
Pp Vo € dT / dx- d( à dT / dx ) / dx = - Qu ©p - Qnoz Ono2 (13) 


with œp the rate of reaction of the propellant, as in Eq. (5), and @no2 that of NO; as in the second term of 
Eq. (12). The summation of Eqs. (5,12,13) through the condensed phase to the surface leads to: 


à dT / dxļ|s = Pp Vb ( Cg Ts - Cp To - Qs ) =Pp Vo Qe (14) 
Qs = Qa + Qnoz (Ynoz,i - Ynoz,s ) (15) 


The first equation is the heat balance at the surface, it means that the heat flux from the flame in the gas 
phase allows the heating and pyrolysis of the propellant. The net heat of decomposition of the propellant 
Q; is exothermic to the extent that some NO), already reacts exothermically in the condensed phase. 


From thermocouple traverses such as that of Fig. 6 and the balance of Eq. (14), the net heat Q, can be 
estimated (again scatter should be expected). The results from various sources are given on Fig. 8. The net 
heat of decomposition is seen to be exothermic and increasing with burning rate (due to an increase in 
pressure). Summation of Eq. (12) yields (with Eq. (5) taken into account) 

Pp Vo ( Ynoai - Yxozs) = Anoz Í (p M/R T ) Ynozexp( - Enoz/ R T ) dx 

Axo: (p M/RT;) Yno: exp( - Enoz/ R Ts) react. 

and with Eq. (8) 


Yno2ii = Yno2s ~ (Pp M/RT,) exp( = Eno2/ R Ts) / exp( = E./RT;) 
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Figure 6: Temperature Profile in the Condensed Phase. 
(Double Base Propellant) 


Example of thermocouple 
Pt — Pt / Rh soldered end to end 


30) Burning rate (mm/s) 


=0.8 10° cm°/s 


+a Microthermocouple ONERA (1100 cal/g) 
A Microthermocouple SUH 
© Reference, microthermocouple KUBOTA 
© + additives, microthermocouple KUBOTA 
e Microthermocouple ZENIN Propellant 
o "Light pipe method" SELEZNEV | "N" 
e "Thermal noise method" (reference) 
m Microthermocouple DENISYUK T,=120°C 


— Corresponds to A,=110" s”, E.=40 10° cal/mole 
Q,=100 cal/g, T,=20°C, d, 


p in atm is indicated 


20;—\ - 4 2 
10+- 
5 
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3 + 
2- coca ee E 
es 
* } 
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Figure 7: Pyrolysis Law. 
(Double Base Propellants) 
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Figure 8: Heat Evolved in the Condensed Phase. 
(Double Base Propellants) 


This relation indicates that the amount of NO, reacting in the condensed phase will increase with pressure, 
and thus Q, will increase, if the reaction rate for NO» catches up with the decreasing residence time in the 
degradation layer, 


Treact. ~ 1 Ivy ~ 1 / exp( = E./RT;) 


(due to Eq. (9)). A plausible law is obtained for a first order NO» reaction with an activation energy Eno2 
of about 15 kcal/mole. 


It should be observed that the heat evolved in the condensed phase Q, has to be affected by the amount of 
NO> present in the propellant, that is by the heat of explosion. 


One important feature of Fig. 8 is that the heat evolved in the condensed phase is not affected by the 
presence of additives. 


4.0 Flame Zone 


As seen above, the reaction between NO, and aldehyde starts in the condensed phase in such a way that 
the surface gas composition [7] indicated on Fig. 5 is obtained, with NO, being significantly less (0.255) 
as compared to what results from the initial degradation of the propellant (0.467, in the case of an 
1100 cal/g propellant) and with NO already present (0.14). The mole balance of Eq. (11) allows to match 
as well the gas analysis at the end of the primary flame ([7], measurements at 11 atm) 


NO: + CH20 — NO + 1/2 CO + 1/2 CO; + 1/2 H20 + 1/2 H? (11) 


Qno2 = 1040 cal/g of NO». 
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It is likely that aldehyde-aldehyde reactions also occur following (again to match the gas analysis results): 
CH,0 + CH20 —> CO + 1/2 CoH4+H20, Qaa. = 389 cal/g of CH20 (16) 
(CHO), + (CHO) — 4CO+2H2, Quataz = - 533 cal/g of (CHO), (17) 


At low pressure when probably only the reaction of (11) can take place, an energy balance between initial 
temperature and end of the primary flame yield: 


Cg Tip - Cp To = Qa + Qnoz Ynoz,i, p < 1 atm (18) 


since the initial NO; is totally consumed in the condensed phase and the primary flame. An evaluation of 
Ty = 1340 K results. It is seen on Fig. 9 that measurements with small thermocouples indicate a large 
increase with pressure of the primary flame temperature from this value. The aldehyde reactions of 
Eqs. (16,17) do not produce energy in significant amount. It is assumed [5] that the NO already present at 
the surface as well as that produced from the NO2-aldehyde reaction react with the layer of carbon residue, 
attached to the surface, which is observed by direct visualization under combustion and after extinction by 
scanning electron microscopy. 


Measurements by microthermocouples 
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Figure 9: Primary Flame Temperature of Double Base Propellants. 


The NO/carbon reaction had been investigated in [11]. If one makes use of the results obtained, 
the conservation of the species carbon can be written as (knowing that the carbon layer regresses with the 
surface) 


Vp dpe / dx = ( Ms / Mro ) ®no (19) 
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vo = - Axor exp( - Exoic/ R T ) pe Sse Yno p M , in g/em’s 
(this form resulting from the way the data of [11] is cast), the reaction balance being assumed to be: 
NO +C = 1/2 N: + CO, Qnoic = 1600 cal/g of NO (20) 


and S,. being the specific surface area of the carbon (at most 10° cm? /g). Reference [11] produces 
(after rounding Eyoic) 


Anoc = 2 10° mole/s cm? atm, Enoc = 30 kcal/mole 


A rough estimate of the amount of N, produced by reaction (20) through the primary flame is given by 
(with m = pp Vp the mass flow rate) 


m Yypfp = = ( 1 /2 ) ( M./ Mno) <No Xf 


<@no> being evaluated at average values through the flame. With p = 11 atm, v, = 0.28 cm/s, 
<Yno> = 0.2, <p> = 0.5 g/em’, <T>= 1400 K and x; = 400 um (from thermocouple measurements) it is 
obtained 


Yno,fp = 0.06 


a reasonable value (with respect to the result of Fig. 5). This tends to indicate that the NO/carbon reaction 
has a kinetics indeed fast enough with respect to the residence time allowed in the primary flame. 


An energy balance taking into account the NO/carbon reaction is written: 
Cg Tf - Cp To = Qa + Qnoz Yno2,i + Yno,cons. Qnorc + Qala Y ala,cons. (21) 


(the consideration of the aldehyde reaction cools the flame by about 80 K). In this balance Yyo,cons. is the 
amount of NO consumed in the primary flame: 


YNOo,cons. = Ynoz,i ( Mno/ Mnoz) - YNO,fp 

in the case of the example of Fig. 5 and taking into account the uncertainty on the measurements: 
YNo,cons. = 9.035 to 0.047 and Eq. (21) results in Ty, = 1420 to 1475 K, an admissible value when compared 
to the results of Fig. 9. 


The temperature profile in the primary flame is controlled by the conservation of energy [5]. 


Mm Cg dT / dx - d(Ag dT / dx) / dx = - Qnoz @Nno2 - Qnoic Ono - Qaid-ald. (22) 


Table 5: Values Considered as Representative for the Gas Phase (Double-Base Propellants) 


cg = 0.35 cal/g K, 
Ag = 1.25 10°( T / 700 )°’ cal/cm s K 


M =30 g/mole 
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The evaluation of the heat flux received at the surface, which will then control the burning rate according 
to Eq. (14), results from the summation of Eq. (22) through the flame zone. Only a true numerical 
evaluation of the temperature and species profiles will give the proper heat flux. However conclusions can 
be drawn from approximate relations. The activation energies of the reactions of Eq. (22) being moderate, 
the flame is distributed and an approximation of the temperature profile is written as: 

(Tr- T)/(Te- Ts) = exp( - 3 x / xe) (23) 


which will produce a shape such as that of Fig. 6. The coefficient of 3 is such that, when x = x;, T;— T is 
5 % of Tt -T,, that is close enough to the final temperature. 


With Eq. (14) 

Ags dT / dx|, = qs = m Qe = Pp Vb (Cg Ts - Cp To - Qs) (14) 
Eq. (23) yields 

qs =3 Ags ( Te- Ts)/ Xf, X#=3 Ags ( Te- Ts)/ MQ. (24) 
The summation of Eq. (22) through the flame results in 

Mm Cg ( Tt- Ts) + qs =3 Qg <g> Xf (<g> positive is the average rate) 

or with Eqs. (14,24), and taking into account an overall equation for conservation of energy 

Cg ( Tr- Ts) = Qg - Qe, 

m=[3<0 Ags ( Te- Ts )/ QJ? (25) 


At very low pressure, ~ 1 atm, when only the NO, reaction probably takes place, the burning rate follows 
pressure according to, see Eq. (12), 


m = Pp Vo ~ (<@no2> )'? ~ p”” exp( - Eno2/ 2 R Tọ) (26) 


a pressure exponent which is indeed observed, see Fig. 2. As the pressure increases the NO/carbon 
reaction takes on more importance and, referring to Eq. (19), 


m = Pp Vp ~ [ p exp( - Eno2/ R Tẹ) + ~ p exp( - Enorc/ R Tẹ) ji (27) 


which, with the increase of the flame temperature with pressure, see Fig. 9, accounts for the pressure 
exponent of 0.7. 


At higher pressures, above about 150 atm for the 1100 cal/g propellant for example, a change in the 
pressure sensitivity, Fig. 2 is observed. This tends to indicate that the secondary flame, probably 
characterized by a second order, with respect to pressure, reaction for NO, comes into the primary flame 
and progressively dominates it, with a pressure exponent, according to Eq. (25), increasing to close to 1. 


Referring to Fig. 10, for the 1250 cal/g propellant, the temperature sensitivity is indicated: 


6, = ( d Inv,/ dT, ) at p given (28) 
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Figure 10: Temperature Sensitivity of a Double Base Propellant. 


According to Eq. (27) the burning rate is under the influence (noting that the condensed phase energy 
balance also contains the effect of Tọ see Eq. (14)) of a premixed flame heat flux and therefore very 
sensitive to changes in the primary flame temperature, in the pressure domain when the two flames are 
separated. From Eq. (21), any change in initial temperature will affect the primary flame temperature and 
therefore induce a change in burning rate. As the pressure rises the primary flame temperature increases 
(up to 1800 K) and then for higher pressures the burning rate comes under the influence of the final flame 
(with a temperature reaching 3110 K for the 1250 cal/g propellant). It is seen from Eq. (27) that a given 
change in To and therefore in T¢ has a smaller impact on the burning rate for higher flame temperatures, 
that is for higher pressures, a tendency observed on Fig. 10. 


5.0 Active Binders 


Various types of active binders, based on nitrocellulose or an inert binder and nitroglycerin or less 
energetic liquid nitrates, can be used, see reference [10] for a complete description. The cross-linked 
double-base binders (XLDB) will be considered here, in which the polymer is cured with an isocyanate 
after mixing with NG. 


Table 6: Values for a XLDB Binder 


Composition: ~ 2/3 NG, 1/3 polyethylene glycol. 
Heat of explosion: 850 cal/g. 
T= 2000 K. pp = 1.42 g/cm’. 
cp = 0.46 cal/g K. 


Ap = 3.9 10* cal/cm s K. dp = 0.6 10° cm’/s 


Ynozi = 0.421. Qa = - 150 cal/g 
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Although the burning rates of the different active binders can be, for a given heat of explosion, somewhat 
different at low pressures [10], above 10 atm the differences become small , see Fig. 11. In the case of a 
double-base propellant and of a XLDB binder (that of the above table), with nearly the same heat of 
explosion, the burning rates for a large range of pressure are very close, Fig. 12. 


o DB ONERA - DYNAMIT NOBEL - HELLER 
e XLDB SNPE 
m XLDB, Ref.10 


an 4 V, (mm/s) (at 100 atm) 


25 
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Heat of explosion (cal/g) 


0 
800 1000 1200 1400 


Figure 11: Burning Rate vs Heat of Explosion. 
(Double Base Propellants and XLDB Binders) 
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XLDB binder 850 cal/g 

+ Lead salt + carbon black 


100 f V, (mm/s) 


OOo 


P (atm) 


1 10 100 1000 


Figure 12: Double Base Propellants and Active Binder Burning Rates. 
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Measurements for XLDB binders of the degradation kinetics, of the surface temperature, of the heat 
evolved in the condensed phase, as in Fig. 8, and of the primary flame temperature, Fig. 9, show that these 
characteristics are very close to those of DB propellants. 


Table 7: Gases Evolved from the Surface, Mass Fractions. XLDB Binder [7] 


Gas analysis at the surface gives results which are qualitatively comparable to those of DB propellants, 
Fig. 5. 


It can then be stated that the combustion mechanisms of the various types of double-base propellants and 
of the active binders are very similar. 


6.0 Mechanisms of Action of Additives 


The incorporation of a few per-cents of lead (and copper) salts and carbon black enables to obtain super- 
rate effects followed by mesa or plateau effects in the burning rate versus pressure laws of double-base 
propellants, Fig. 3, as well as of active binders, Fig. 12, although in the latter case these effects are much 
less pronounced. 


What is thought to be the mechanism of action of the additives has been presented in references [5,6] 
by the authors and in reference [12]. 


It has been found that the active part of the lead salt is the oxide of lead which accumulates above the 
propellant surface, after the salt has been trapped in the carbon residue layer which can be observed, 
immersed in the primary flame (the decomposition kinetics of the salt is slower than that of the propellant 
components and it thus emerges from the surface unchanged). If the propellant (when its heat of explosion 
increases) or the active binder naturally produces less carbon residue, then the lead salt particles are in 
large part ejected from the surface and cannot act. This is the case probably when the amount of 
nitrocellulose is reduced (hot double-base propellants) or almost absent (XLDB binders). In the latter case 
however the inert polymer leaves some carbonaceous residue. Addition of carbon black is probably 
favorable because it accumulates on the surface in the naturally produced carbon layer. 


It has been found by the authors that PbO reacts preferentially with aldehydes to form carbon and CO). 
It has been observed systematically [6] that there is a relation between the amount of carbon residue and 
the importance of the super-rate observed (for example, depending on the fabrication process: solventless 
extrusion, powder casting, mixing). It was seen that the primary flame is due to an NO,-aldehyde reaction. 
But NO starts to react with carbon as well close to the surface, leading to the increase in primary flame 
temperature of Fig. 9. It is believed that the extra carbon produced in the presence of additives enhances 
the NO/C reaction (see ref. [12] for more results on this reaction), depositing extra energy in the primary 
flame (an increase in primary flame temperature is observed in the presence of additives [5]) and resulting 
in a higher heat flux to the surface, and thus a higher burning rate. 


Figure 13 shows surface structures of a cool propellant with a strong super-rate and a corresponding thick 
carbon residue. In this case, visualization shows that a physical effect occurs in which the secondary 
luminous flame attaches in streaks to the carbon layer and deposits its high temperature (~ 2200 K as 
compared to the ~ 1400 K of the primary flame) closer to the surface. The abrupt end of the super-rate, 
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the mesa effect, is believed to be due to the too thick carbon residue being expelled from the surface. 
At the end of the mesa effect the surface is almost clean of carbon. 
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Figure 13: Super Rate Effects in Double Base Propellants. 
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In the case of hot propellants, Fig. 3, the super-rate is probably due only to the chemical effect of the 
enhanced NO/C reaction. As the pressure increases the secondary flame, where NO will react anyway, 
merges into the primary flame and progressively the modified propellant is caught up by the reference 
propellant, a plateau effect thus resulting. 


In the high pressure domain when the flame system has reduced to one overall flame, a second super-rate 
occurs (when only lead salt or oxide is added), see Fig. 13 for the cool propellant and Fig. 3 for a 
1250 cal/g propellant. This effect seems to be purely physical, related to the presence of lead oxide 
particles accumulating and imbedding into the surface, Fig. 13, with for example an enhancement of the 
thermal conductivity of the flame zone or more probably a flame holding effect (the protruding particles 
perturb the flow from the surface). At higher pressures and burning rates the thickness of the condensed 
phase heated zone and reaction layer and flame zone becomes so small that the particles will not attach to 
the surface or will be too large to perturb the combustion process. 


PYROLYSIS OF INERT BINDERS 


A number of books and works has been devoted to the behavior of polymers, whether or not usable as 
binders, under thermal loads, references [13-17] are examples. Much work has been carried out with 
thermogravimetric analysis (TGA) or differential scanning calorimetry (DSC), with heating rates at most 
of the order of 1°C/s. Under linear pyrolysis (for a binder within a solid propellant) the rate of temperature 
increase is of the order of 10° °C/s. It is far from obvious a priori that the degradation kinetics will remain 
the same. In Ref. [8] it was attempted to establish that this is indeed the case for a number of polymers. 


Although it is hardly a propellant binder, Teflon is an interesting reference polymer. Its degradation 
kinetics (obtained by TGA) and thermal properties [8] ( Ap = 6.34 10% cal/em K s, pp = 2.1 g/cm’, 
Cp = 0.25 cal/g K) are indicated on Fig. 14. In order to extrapolate these characteristics to the regime of 
linear pyrolysis (obtained experimentally by pressing the sample on a hot plate) the procedure of 
reference [8], also explained in the condensed phase paragraph of the double-base propellants chapter, 
is applied. In the case of a first order (with respect to the non degraded polymer) reaction, the relation 
between regression rate and surface temperature is (again, numerical computation shows this relation to be 
accurate to about 1 %) 


vi = (dp / &) Ac exp(- če ) /[ ( - INYps) (1 - To/ Ts - Qs/ Cp Ts) +Qs/ ep Ts ] (1) 
& =E./RT; 

In this relation Q, is the heat evolved in the condensed phase, in this case endothermic and equal to 
- 340 cal/g ( Yp, mass fraction of the remaining polymer at the surface can be set at 0.01 ). It is seen on 
Fig. 14 that there is a good match between extrapolated law and measurements. These measurements are 


obtained under various atmospheres, showing no influence of this factor. The conclusion is then reached 
that the pyrolysis of such a polymer is an irreversible thermal mechanism. 
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Figure 14: Linear Pyrolysis and Extrapolation of Thermogravimetric Results (Example of Teflon). 


In the case of an actual, widely used, propellant binder such as hydroxyl terminated polybutadiene (HTPB) 
the same extrapolation can be made, Fig. 15 (the thermal properties used: 4, = 3.6 10“ cal/cm K s, 
Pp = 0.92 g/cm’, Cp = 0.39 cal/g K), and compared to the results of Ref. [15], obtained in a hybrid motor, 
gaseous oxygen upon HTPB, with the surface temperature being measured by infra-red pyrometry 
(with some dispersion). A measurement obtained by the authors with CO, laser heating and IR pyrometer 
is also indicated. The agreement, considering the extent of the extrapolation, is very satisfactory. 


Further on, the pyrolysis law based on the TGA and DSC kinetics will be considered to apply. 


RTO-EN-023 


Combustion of Solid Propellants 


— From degradation kinetics 
o IR measurements by Strand, Ref.15 
e CO, laser heating, 
I.R. pyrometer, ONERA 


Mass flux (g/cm’s) 
a | 


E, = 48.6 kcal/mode 
(T. Brill = 46) 

—| A, =2.63 10" s" 

d, = 110° cm’/s 

Q; = -450 cal/g 


0.10 = 2 


0.01 
7.00E-4 9.00E-4 1.10E-3 
Surface temperature 1/T (K) 


Figure 15: Pyrolysis Law for HTPB. 


Equation (1) and Fig. 15 indicate how the surface temperature adjusts itself to allow the polymer to 
degrade into gases when the regression rate changes. Further considerations to better know the binder 
behavior are the nature of the gases resulting from the pyrolysis and the corresponding heat of 
degradation. 


The heat of ablation of HTPB was measured in Ref. [16] by relating the mass ablated to the radiation heat 
flux received by a sample. For a regression rate of 0.4 mm/s (and from Fig. 15 T, ~ 1060 K) it is obtained: 


havi = Cp ( T; = To) = Q; =750 cal/g 
Q, = -450 cal/g 
It is found that this value is compatible with the production of mostly C4Hį as pyrolysis gas. Gas analysis 
at low temperature [13] reveals a complex set of gases (~ 45 % in mass butadiene, ethylene, propene...) 


resulting from the degradation of HTPB. If one evaluates the difference in heats of formation between the 
initial material and butadiene it is found: 


Ah = h° (C, H, ) - h° (HTPB) = 496 cal/g - 5 cal/g = 490 cal/g 


close to the measured heat of degradation. At high heating rates butadiene is probably the major 
degradation gas. 


COMBUSTION OF AMMONIUM PERCHLORATE 


Ammonium perchlorate (AP), NH; C10;,, is a widely used oxidizer and as such has been the object of 
numerous investigations. References [18 to 22] are a sampling, with [22] giving a detailed list. A view of 
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the combustion of AP is presented herein attempting to make use as much as possible of the various 
experimental data available. Due to the large number of works on AP combustion, somewhat 
contradictory interpretations and corresponding models have been produced. A simplified model is 
presented, which is considered to represent reasonably the combustion mechanism of AP, although it will 
not be in agreement with all of the above mentioned interpretations. 


1.0 Condensed Phase Behavior 


The considerations presented previously for the condensed phase of a pyrolysing monopropellant apply to 
AP. The conduction zone has a thickness 


cond. = ( dp Vp ) In10° 


(with the thermal diffusivity [21] dp = 1.2 10° cm’/s at an average temperature in the heat wave), 
thus equal to = 55 um, for a burning rate of 10 mm/s. 


Table 8: Condensed Phase Values for AP [19,21] 


Cp = 0.31 cal/g K (orthorhombic phase < 513 K) = 0.365 cal/g K (cubic phase) 
Pp = 1.95 g/em’, dp = 2.5 10° - 4.55 10° T(°C)cm?/s 


Based on observations by scanning electron microscopy after extinction, the idea has been advanced that 
the self-deflagration of AP, possible only above 20 atm, requires that a large amount of exothermic 
reaction already takes place in the condensed phase in a thin liquid layer (above a melting temperature 
estimated at 835 K) [19,20,21,22]. 


2.0 Energetics of the AP Combustion 


The model of Ref. [19] is subscribed to in order to describe the combustion of AP alone. The AP 
undergoes a phase transition at 513 K, melts around 830 K and, in the thin (a few microns) superficial 
liquid layer thus created, an exothermic reaction, affecting 70 % of the AP, takes place and creates the 
final combustion gases, O, in particular. The remaining 30 % of the AP sublime into NH; and HCIO, 
which react exothermically in a premixed flame very close to the surface (a few microns), Fig. 16. 


1.015 O,+.265 N,+.12 N,O Final composition 
+.23 NO+1.62 H,O+.76 HCI+.12 CI, (After FRIEDMAN) 


i 


AP premixed flame (30%) 
Q, = 265 cal/g 


Condensed phase AP 
reaction (70%) 
Q, = 120 cal/g 


Figure 16: Autonomous Combustion of Ammonium Perchlorate. 
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From the data collected in [19] the change of enthalpy per gram of AP required to heat up the AP to its 
surface temperature, Ts ap, is estimated as 


A hyap = 0.31 ( 513 - 293 ) + 21 + 0.365 ( 835 - 513 ) + 60 + 0.328 ( Ts. ap - 830 ), 
= 266 + 0.328 ( Ts. ap - 835 ) cal/g (1) 


where the heat of transition and the heat of liquefaction appear. The exothermic condensed phase 
degradation involves an enthalpy change of [19] 


A hp ap ~ - 380 cal/gram of reacting AP. (2) 
The heat of sublimation is 58 V 2 kcal/mole or 
Ah s ap = 476 to 510 cal/gram of subliming AP. (3) 


The adiabatic flame temperature for the combustion of AP alone has been estimated in [18] to be 
T; ap’ = 1205 K, corresponding to the reaction 


NH; C10, — 0.265 N2 + 0.12 N20 + 0.23 NO+ 1.62 H2O + 0.76 HCl + 0.12 Cl, + 1.015 O2 (4) 
The change of enthalpy, Ahc.ap, corresponding to the combustion of the sublimed NH; and HCIO, (to give 
the combustion products of (4)) is obtained from the equation expressing the conservation of enthalpy 
between the unreacted AP at initial temperature and the combustion products downstream of the flame: 

0.3 Ahcap + Cy ( Tiap“ - Tsar) + 0.3 Ahsap + 0.7 Altp,ap + Ahy ap = 0 (5) 
For cg ~ 0.3, value taken in [19], it is found 


Ahc ap = - 850 to - 885 cal/gram of reacting AP, (6) 


depending on the value adopted for the heat of sublimation, 476 or 510 cal/g, and independently of the 
value of Ts,ap in the range found in [19]. 


It should be noticed that the value given here for the transformation of the AP into gases (that is NH; and 
HCI1O, for the subliming 30 % and the combustion gases of (4) for the 70 % reacting in the condensed 
phase), namely 


- Qs = 0.3 (476 to 510 ) + 0.7 ( - 380 ) = - 123 to - 113 cal/g of AP, (7) 


is also found in [20], where a model for the combustion of AP similar to that of [19] is adopted. The heat 
evolved in the condensed phase Q, ( > 0 if exothermic) will be set equal to 120 cal/g. 


One last check of consistency can be performed: with the enthalpy of formation of AP at h’,p = - 602 cal/g 
and that of the combustion products of (4) h°cp = - 877 cal/g, an overall energy balance between initial AP 
and combustion products is: 


d, 10 0 
Cg Trap + bcp = Cortho To +h ap (8) 


resulting in T;ap “’ = 1215 K, close enough to the previous value. 
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3.0 Surface Pyrolysis of AP 


Attempts to measure the surface temperature give values between 670 and 973 K. These results, obtained 
either by using thermocouples imbedded in the AP pellet, or by measuring the radiation emitted by the 
surface, are always associated with some uncertainty due to the operating methods. In effect, the large size 
of the thermocouples, in relation to the temperature gradients encountered, favors errors; further, 
the measurement represents an averaging of surrounding conditions. The temperatures thus obtained are 
therefore probably somewhat lower than in reality. In the case of experiments using an optical technique, 
the temperature is deduced from measuring the radiation emitted by the surface and transmitted through 
the gaseous layer. The disturbance caused by the radiation of the gases and the screening action of the 
flame limits the application of this method to 60 atm. That is why the original technique suggested by 
Seleznev [23], carried out by using a sapphire light guide inserted in the solid substance and by reading the 
infra-red emission in the direction of the condensed phase, has the considerable advantage of providing a 
direct measurement of the radiation emitted by the surface, without any hot gases and the reaction of the 
flame being interposed. Its application can therefore be extended to high pressures and the measurements 
appear to be more convincing. These results enable to determine the AP pyrolysis law. 


Further, estimates of the AP melting temperature are mentioned by a number of authors; the values 
suggested vary from 715 to 865 + 20 K and are useful in interpreting the extinction phenomenon at low 
pressure. 


The pyrolysis law of the AP is written as: 


Map = Pap Vo,ar = Asap €XP( - Es,ap/ R Tsar ) (9) 


The parameters characterizing the pyrolysis law are determined so as to obtain a good agreement between 
the rates and the surface temperatures measured by Seleznev [23] (Fig. 17). The activation energy 
obtained is 20 kcal/mole, a figure compatible with the various estimates encountered. The measurement of 
the surface temperature made at 40 atm is the only one which deviates from that computed by the 
pyrolysis law used. On the other hand, for the critical rate of 0.27 cm/s obtained at 20 atm (AP combustion 
pressure limit), this law allows for a surface temperature of 830 K, corresponding to the assumed AP 
melting temperature [19]. 


Es 4p = 20 kcal/mole 
Asap = 96000 g/cm’-s 
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Figure 17: Pyrolysis Law for Ammonium Perchlorate. 
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4.0 Flame Structure of the AP Combustion 


The approach of [19] considers 14 different reactions to describe the flame zone of the AP, involving 30 % 
of the material sublimed into NH; and HC1O,. An overall second order reaction is then obtained 


NH; + HClO, > products of Eq. (4) 
with an activation energy of E, ap = 15 kcal/mole. 


A simplified model is constructed like the one used in the chapter on the double base flame. 
The temperature profile is approximated as 


(Trap - T)/( Tear“ - Ts,ap ) = exp( - 3 x/ xr) (10) 
with the boundary condition 

qs = Ags dT / dx|s = map Qe (11) 
Q. = Ahyar - Qs 

Q, is the heat required to bring the AP from To to gases at the surface, see Eqs. (1) and (7). 

Combining (10) and (11) results in 

Map = 3 Ags ( Tear” - Tsar) / X£ Qe (12) 


Summation of the energy equation from T;.ap, x = 0, to x¢ yields, after using Eqs. (11) and (12) and Eq. (5) 
written as 


Qg + Cy ( Trav - Ts ap) +Q. = 0 

map = [ 3 <@> Ags (Trap - Tsar ) / Qe ]'” (13) 
In this equation the average reaction rate is written as 

3 <@.> = p° Agar exp( - Egap/R Trap“) (14) 


expressing a second order overall rate. The burning rate is then found to follow pressure with an exponent 
close to 1. Table 9 indicates various values for AP combustion. 


Table 9: Values Considered as Representative for the AP Flame Zone 


Pap = 1.95 g/cm? 

Esap = 20 kcal/mole 
Asap = 96000 g/cm” s 
Egap = 15 kcal/mole 
Agap = 700 g/cm’ s atm? 


y =30% 
Cg =0.3 cal/g K 
À = 1.9 10* cal/cm sK 


g 
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Taking into account the various values given, it is obtained for vp = 10 mm/s ( Ts.ap ~ 925 K ) the flame 
stand-off distance, Eq. (12), ( Q. = 175 cal/g ), x= 5 um, the height to reach Tp within 5 %. 


The only input data which is not defined from outside considerations is the pre-exponential factor Ag ap 
used as a floating parameter for the model. Various numerical values associated with the input parameters 
of the model are brought together in the above table. Just by the choice of the prefactor A, ap adjusted at 
700 g/cm* s atm’, the model satisfactorily reproduces variations in the AP burning rate due to pressure as 
well as due to the change in the initial temperature of the product. Figure 18 provides a comparison of 
computed rate curves with experimental points [21]. 


| V, critical = 2.65 mm/s | 


| V, (mm/s) 7 A V, (mm/s) | 
a +: 
10t al A 10} “ 4 
7 [Pn= 202 bar) — | [Pa on = 10.1 bar | 
j | 
od * 298 K + 323 K 
1 _ 1 , diiit 
z REEE: 
10} ! 10 Pa J 
Poi = 11.9 bar | f Z Pa = 14.6 bar 
| "4 b 4 
< 423 K - 373K | 
4| P (bar) r P (bar) | 
10 100 10 100 


(Exp. at 25, 50, 150°C from T. Boggs) 


Figure 18: Burning Rate of Ammonium Perchlorate at Different Initial Temperatures. 


At sufficiently high pressures, the energy transmitted from the flame toward the surface, to which must be 
added the effect of the superficial exothermic reactions, maintains the surface temperature above the AP 
melting point. When the pressure falls, the premixed flame moves away and the surface temperature can 
then drop below the limiting value, thus causing the disappearance of the liquid surface layer which was 
enabling the exothermic reactions to occur. The energy from the flame is then much too small to maintain 
a pyrolysis which has become strongly endothermic and AP no longer burns. This minimum pressure, 
beyond which the combustion cannot propagate itself, sets the pressure limit for AP self-degradation. 


On the basis of this hypothesis, the combustion pressure limit is reached when the surface temperature is 
equal to the AP melting temperature. It is interesting to use the model in order to follow the variation of 
the computed limiting pressure with the initial temperature. At atmospheric temperature, the combustion 
limit is 20 atm, in agreement with experimental results. The critical rate of 0.27 cm/s and the surface 
temperature of 830 K, representing the AP melting temperature, corresponds to it. The computation 
method consists, for the initial temperature varying between 0 and 200°C, in finding for what pressure the 
surface temperature is equal to 830 K. A comparison between the computed and experimental pressure 
limits is very good, Fig. 19, and confirms the soundness of the hypothesis following which AP only burns 
if the surface temperature exceeds its melting temperature. 
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Figure 19: Limiting Pressure for Combustion of Ammonium Perchlorate. 


One further set of results is presented. Questions have been raised about the combustion of AP at high 
pressures, above about 100 atm, with many conflicting results. Measurements of burning rates on carefully 
inhibited samples are presented on Fig. 20. It is concluded that no strong change of combustion regime is 
observed. The model based on the above presented hypothesis seem to follow quite well the results to high 
pressures. 


o Inhibited samples (ONERA, 1988) 


— Computed burning rate 
Burning rate (mm/s) 
100 
50 F 
10 [ 


Pressure (atm) 


1 
10 50 100 500 1000 


Figure 20: Burning Rate of Ammonium Perchlorate at High Pressures. 
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COMBUSTION OF HMX 
The combustion of HMX, which has balanced 
oxidizing and combustible elements (to reach CO, 0, p 
H:O and N3) and a combustion temperature of 3280 N 
K, is controlled by processes in the condensed CH,—N 
phase as well as in the flame zone. Information Ou ai? 
about these processes is to be found in references "o 7; 
[24 to 31] (also [21)). HC 

Ny- cHy 

N, 
J `o 


Figure 21: Chemical Structure of HMX. 


1.0 Condensed Phase Processes 


The kinetics of the decomposition of HMX can be obtained by differential thermal analysis [24,29]. 
A fairly extensive review is given in Ref. [25]. Results obtained around the (assumed) melting temperature 
of 556 K are indicated on Fig. 22, with the kinetics extracted for a reaction order of 1 (refer to [29] for the 
method utilized). The condensed phase properties considered to be representative are indicated below. 


30 mg, 1 atm 
Order n = 1 
Activation energy 
E, = 50 kcal/mole 


A, = 0.7 to 1.6 10” s” 
10 t Heating rate (°C/mn) 


dQ/dt 


1 bee 
0.1 ee E 
1.81 1.86 1.91 1.96 
(552 K) (523 K) 10°/T (K) 
Decomposition peak 
temperature 


Figure 22: Differential Thermal Analysis of HMX Decomposition. 


RTO-EN-023 


NATO 


OTAN 


Combustion of Solid Propellants 


Table 10: Condensed Phase Properties (Average in the Thermal Wave) for HMX [31] 


Pp = 1.9 g/cm’, cp = 0.33 cal/g K 


dp / Àp / Pp Cp = 1 10° cm/s 
(also found from thermocouple traverses) 


Ap = 6.3 10% cal/em s K 


Another technique to obtain the decomposition kinetics, at somewhat higher temperatures, is through 
ignition tests, by exposing the sample to a given surface heat flux and detecting the delay for the first 
exothermic ignition reaction (rapid deviation of the surface temperature from that of an inert material). 
The results are seen on Fig. 23 and confirm those obtained by DTA (the first method gives a good estimate 
of the activation energy, whereas the ignition experiments are better for estimating the prefactor). All these 
results correspond to an irreversible thermal decomposition of the HMX bonds with no influence of the 
pressure. 


| Correlation 


| E, = 50 kcal/mole 
| A, = 1.45 10” s" 


T = (A, P, C)” = 0.02 cal/K cm’ s” 
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Figure 23: Ignition of HMX. 


Finally, the investigation of the condensed phase processes under combustion is carried out by 
determining the temperature profile through the combustion wave by use of micro-thermocouples. As was 
seen on several occasions the thermal wave thickness is: 


Ccond, = In107 dy / Vp (1) 


that is, in the case of Fig. 24, econg, ~ 220 um. The thermocouple junction has to be very small, in the 
present case 5 um platinum wires are welded end to end by electric discharge with the junction at about 
this size. A measurement of the thermal diffusivity is also obtained, close to that indicated in Table 10. 
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The relationship between the surface temperature (obtained from measurements such as that of Fig. 24) 
and the burning rate is displayed on Fig. 25. Scatter of the measurements is hard to avoid (about 50 K). 
Also indicated are thermocouple results from reference [28]. The pyrolysis law is also established from the 
decomposition kinetics obtained by DTA and by ignition experiments, making use of the approach 
mentioned for double-base propellants as well as for inert binders, according to which [8], for a first order 
reaction, 


Ve = ( dp/ c) Ac exp( - & )/[(-InYps) (1 - To/Ts- Qs/cpTs) + Qs/ cp Ts] (2) 


& = E. / R T, and Q; is the heat evolved in the superficial reaction layer of the HMX; its value will be seen 
next. The amount of HMX at the surface Y, is set at 0.01. It is seen that there is continuity between the 
decomposition (thermal breaking of the chemical bonds, with probably no participation of the vaporization 
of the HMX) under DTA conditions, = 550 K, for ignition, from 600 to 700 K, and under combustion, 
up to 900 K. This conclusion on such a continuity, which is not a priori guaranteed, was also reached for 
double-base propellants (as well as for various binders). 


i 32 e 7; measurements 
V, = 2.1 mm/s — d, = 1.110°cm*/s o Mitani, Ref.28 


HMX + wax (2%) —— Pyrolysis law 


p = 11 atm hibura rars aiid E, = 50 kcal/mode 
m/s R 
10! g 9 | ge tas 105" 


Temperature platea F | | =110° 
4 T(K) p paeau | d, = 110° cm’/s 
2000 — ) „m 5 
L | 4 }__| ra | o 
1600 | 1 | È 2 — 
ail 1 $ ] r 
1200 —— | a t 
S | | E | | 
a a] Ei 0.5} 
800 4+ 5+, _,__ ++ i 
| all 
f a = T 
ga” T; = 760 K 0.2 oo 
400 fargar 59 um —— 4 
| i | R N ee i Surface temperature 10°/T, 
glll- i Timeinms |- i 1.1 1.2 1.3 14 15 
0 20 40 60 80 100 10K 770K 667 K 
Figure 24: Temperature Profile, Under Figure 25: Pyrolysis Law for HMX. 


Combustion, in HMX by Thermocouple. 


From the temperature profiles, as on Fig. 24, the heat evolved in the condensed phase can be evaluated by 
use of the relation: 


Ag dT / dxļs = Pp Vo [ Cg Ts - Cp To - Qs | (3) 


with Q, > 0 if exothermic. The results of Fig. 26 show that indeed the transformation of HMX into gases is 
exothermic. The energy absorbed by the breaking of the HMX bonds is more than compensated by 
exothermic reactions taking place in the superficial degradation layer, probably between NO, and HCHO 
(see further for the identification of these gases). This mechanism was already observed for double-base 
propellants. 
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e Thermocouple measurements 
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Figure 26: Heat Evolved in the Condensed Phase of HMX. 


Experiments have been performed [7] by maintaining the linear regression (at about 1 mm/s) of HMX 
samples with an external heat flux (radiation or contact with a heated plate), under vacuum so as to avoid 
any gas flame. The sampled gases are then analyzed by mass spectrometry. The following table indicates 
the results. 


Table 11: Gases Evolved from the Surface of HMX (Mass Fractions [7]) 


22.2 17.5 19.1 3 1.4 1.6 3.7 0.1 14.1 16.8 


It is seen that the initial degradation produces probably NO, and N20 in similar amounts and HCHO and 
HCN. See also for such conclusions reference [32]. Exothermic reactions involving NO» occur in the 
superficial degradation layer to give a large amount of NO. Results in rough agreement with those of the 
above table have also been obtained in Ref. [33] with an infra-red analysis technique. 


2.0 Gas Phase Behavior 


The production from the condensed phase of HMX of several oxidizing gases, NO», N.O and NO, 
can create a two-stage flame. This was seen to be the case for double-base propellants for which NO, and 
NO are created in the condensed phase: the primary flame involves NO» and the secondary flame NO; 
above ~ 200 atm the two flames merge into one. In the case of HMX, observation at pressures around 
1 atm reveals the existence of a dark induction zone and a detached luminous flame, similar to those of 
double-base propellants. Also the micro-thermocouple traverses below about 10 atm, such as in Fig. 24, 
show a plateau in the temperature profile at about 2000 K, much below the final temperature of 3280 K. 


The burning rate of HMX, as single crystals or as pressed samples, obtained in Ref. [25] is shown on 
Fig. 27. The evolution of this burning rate with pressure shows that around 20 atm a pressure exponent of 
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about | is attained, revealing that the staged flame has collapsed into one (as in the case of double-base 
propellants above about 200 atm). 


Results of Boggs (NWC) 
| o Pressed samples 20° C 
e Crystals 20°C 
Op pgp, = 0.15 to 0.2% /°C 


P exp. 


A 
Burning rate (mm/s 
100 + g ( Is) 
50 | | 
2 nd order NO,/N,O/NO reaction 
E,, = 25 kcal/mole ee 
n= 
10 | 
5 | 
1st order NO, reaction 
E, = 15 kcal/mole 
n= 0.7 
1 = aes on. ere 
1 5 10 50 100 500 1000 


Figure 27: Burning Rate of HMX. 


It was seen in the chapter on double-base propellants that in the case of a distributed flame an 
approximation for the temperature profile is: 


(Te-T)/(Te- Ts) = exp(-3x/x¢), Xt=3 À as (Tr-Ts)/ MQ., Q. = cg Ts - cp To - Qs (4) 


(at v, = 10 mm/s this will give xs = 48 um, the height to reach Ts within 5 %) and the burning rate 
becomes: 


M = Pp Vo = [ 3 <g> Ags ( Te- Ts )/ Qe T? (5) 
With <œ, > the summed reaction rate through the flame zone. More complete descriptions of the flame 


zone can be found in references [29,30], the conclusions of which are essentially those which can be 
extracted from the above simplified approach. 


Table 12: Values Considered to be Representative of the Gas Phase of HMX 


cg = 0.35 cal/g K, Ag = 1.25 10° ( T / 700 )°’ cal/em s K 


At low pressures the primary flame, involving NO, in a first order reaction, controls the burning rate, with, 
according to Eq. (5), a pressure exponent around 0.5. At higher pressures the collapsed flame is probably 
dominated by a second order reaction involving NO and N20, with a pressure exponent close to 1. 


4-34 RTO-EN-023 


Combustion of Solid Propellants 


COMBUSTION OF RDX 
RDX has a structure similar to that of HMX. The present 
authors have not considered this component very much, NO, 
as opposed to HMX. The information comes mostly from Ni 
the work of Ref. [34]. SZN 

Ha Çh 

N N 
N 
Rf \ Cc’ NO, 
s 2 


Figure 28: Chemical Structure of RDX. 


It has been seen that HMX has a behavior somewhat similar to that of a DB propellant, due to the fact that 
it creates gases which can produce primary and secondary flames, see Table 11. In this case, a very 
simplified model, based on overall kinetics, shows that the secondary flame dominates the burning rate as 
soon as 20 to 30 atm (for DB propellants this occurs above 100 atm), in such a way that there is no 
possibility of creating a super-rate effect (and the accompanying plateau). A much more elaborate model, 
see Fig. 30, has been presented by Yang et al [34] for RDX. This model takes into account the two 
condensed phase decomposition paths suggested by Brill, but also incorporates a large amount of 
evaporation of the RDX (it has to be noted that in the case of HMX, referring to Figs. 22 to 25, 
the conclusion is reached that the condensed phase process under combustion corresponds to complete 
decomposition, as observed under DSC-DTA conditions, and thus with no contribution of evaporation). 
In the melt layer of the RDX (melting temperature supposed at about 470 K) the products of the first 
decomposition start reacting, much as was seen for double base propellants and for HMX, this being 
manifested by the gas composition obtained at the surface, see Fig. 30. At the surface, with a temperature 
depicted on Fig. 29, an important void fraction is computed, corresponding to the amount of RDX which 
has decomposed in the condensed phase. 
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Figure 29: Conditions at the Surface of RDX [34]. 
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From T. Brill 
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Figure 30: Modeling of RDX Combustion [34]. 
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A new caged nitramine called either CL20 or HNIW 
(hexaaza hexanitro isowurtzitane) has appeared in the 
last few years (1987). Because of its high energy 
content it is envisioned that when combined to an 
energetic binder, as compared to HMX or RDX 
compositions, it will be better performing in terms of 
impulsion and level of burning rate, with smokeless 
properties maintained and similar safety behavior, 
Ref. [35]. 
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@ N-NO> 


Figure 31: Chemical Structure of HNIW. 


Differential scanning calorimetry (DSC) and thermogravimetry (TG) contribute to the thermal knowledge 
of the components to get specific heat, degradation steps and kinetic constants of the material. 
Helium pycnometry and a thermal diffusivity set-up, at ambient temperature and atmospheric pressure, 
give the other physical properties. Table 13 brings together the results obtained by these methods. 


See Ref. [36] for full details. 


Table 13: Thermal and Physical Properties of HNIW 


Properties Symbols SI units HNIW 
Specific heat p [J.kg"K"] 1004 (+ 13 %) 
Specific mass (density) p [kg.m™] 1980 (+ 1 %) 
Thermal diffusivity din [m.s] 0.457*107 (+ 6 %) 
Thermal conductivity | A=p*cp*dm | [J.m Ks] 0.068 
Activation energy Ea [J.mole™'] 224000 (+ 4.5 %) 
Frequency factor ko or A [s’] Ifn*=1 > 5*10°” 


*n is the reaction order. 


The gases evolved from the condensed phase are obtained, as was done for the double base propellants 
and HMxX, with the linear pyrolysis technique under vacuum so as to avoid the gas flame, the regression 


rate is about 0.5 mm/s. 
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Figure 32: Set-Up for the Linear Pyrolysis. 
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Table 14 shows the results. NO, and NO are present in large amount, as for HMX, N20 is in much smaller 
amount. This gas composition should lead in principle to a staged flame, but due to the high temperature 
reached at the end of the primary flame, which involves a high mass fraction of NOs, it is likely that the 
two flames will collapse into one even faster than for HMX. 


Table 14: Molar and Mass Fractions of the Products from HNIW Linear Pyrolysis 


Species 


NO, 


CO, 


N20 


HNCO 


HCHO 


NO 


HCN 


H,0 co N2 


Solid Solid 
residue 1 residue 2 


Molar 
fraction 


0.017 


0.024 


0.010 


0.017 


0.142 


0.046 


0.049 0.044 0.029 


Mass 
fraction 


0.018 


0.027 


0.011 


0.013 


0.107 


0.031 


0.022 0.031 0.020 


Mass 
fraction 
Gas 
+Res. 


0.015 


0.023 


0.009 


0.0109 


0.0899 


0.026 


0.0185 | 0.0259 0.017 


Figure 33 gives burning rate measurements for HNIW, on a “strand burner”, but it has been now given up 
in favor of a closed autopressurizable bomb with an ultrasonic measurement technique. This has the 
advantage of giving the evolution of the burning rate versus pressure with a very low number of 
experiments, by excursion of a pressure shift during the experiment. Results on ADN, ammonium 
dinitramide, NH,-N-(NO>)2, another new oxidizer of interest, is also indicated. These results are in good 


agreement with most of the literature. 
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Figure 33: Burning Rate of HNIW (and ADN). 
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Small thermocouples imbedded in the burning sample allow to measure the temperature profile, Fig. 34, 
the surface temperature, in the vicinity of 800 K, Fig. 36, and from the heat balance at the surface 
an evaluation of the heat evolved in the condensed phase is obtained, Fig. 35, it is exothermic around 
400 kJ / kg, about 100 cal/g similar to the value observed for HMX. 
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Figure 34: Temperature Profile in Burning HNIW. 


700000 


600000 


500000 


400000 


300000 


200000 


100000 


0 


Qs(y) evolution versus pressure(x) 
y = 450513x + 140533 
RÊ = 0,8592 


Pressure [MPa]. 


1,2 


RTO-EN-023 


Figure 35: Energy Released in the Condensed Phase of HNIW, J / kg. 
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Figure 36: Surface Temperature of HNIW. 


measured surface temperatures and the correlation obtained from the extrapolation 


(according to the method exposed for HMX) of the condensed phase decomposition kinetics. Taking into 
account a large uncertainty related to the size of the thermocouple, the agreement is appropriate and found 


to mean that the dec 


omposition of the HNIW into gases is a purely thermal phenomenon, as measured by 


DSC, not affected by pressure and with no evaporation entering in the process. 
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Figure 37: Surface Temperature vs Regression Rate Correlation. 
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COMBUSTION OF COMPOSITE PROPELLANTS 


This chapter is devoted to the description of propellants made of an oxidizer, ammonium perchlorate or 
HMX, RDX, HNIW, and a binder, inert such as polybutadiene or active, a mixture of a liquid nitrate and a 
polymer. For the sake of clarity only these two categories will be considered, the main point being, as was 
said in the introduction, to attempt to improve the understanding of the combustion of propellants rather 
than to present an extensive catalog of results. 


1.0 Comparative Picture of Composite Propellants Combustion 


Drawing from the results presented previously for the various components, the comparative picture of 
AP-inert binder and HMX or another nitramine -active binder is found on Fig. 38. 
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Figure 38: Comparative Picture of AP-Inert Binder and HMX-Active Binder. 
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The combustion of AP_results in a premixed flame at about 5 um from the surface at ~ 1300 K or higher. 
From this flame 1 mole of O, comes out for each initial AP mole. From the binder surface at nearly the 
same temperature, combustible gases are ejected, which, after diffusion, react with O2. As will be seen 
further on, the diffusion flame height is directly related to the AP particle size. The smaller the particle 
size the closer is the flame and the higher the heat flux to the surface and therefore the burning rate of the 
propellant. The pyrolysis of the inert binder is purely endothermic (heat required to bring the temperature 
to 1300 K and heat to decompose it into gases). In the case of particles of a few tens of um, the burning of 
the AP being close to adiabatic (its flame receives some heat flux from the final flame, but computation 
results show that there is only a moderate deviation from adiabatic conditions), the heat flux from the final 
flame serves primarily to keep the binder regressing. 


The combustion of HMX or another nitramine is also through a premixed flame, about 50 um from the 
surface, reaching the final stage of 3280 K (for HMX). The gases emitted from this flame cannot sustain 
any further combustion. The active binder goes through its own combustion, with a primary flame 
reaching ~ 1550 K some 50 um from the surface. The final flame somewhat further away reaches about 
2000K. There is no direct interaction between the two components of the propellant. The burning of the 
propellant is then an average of the individual burning rates. There exist however an indirect interaction of 
the active binder on the HMX particles. As was seen, the thermal properties of HMX and the active binder 
are close and the HMX particles are immersed in the temperature profile of the binder. Upon reaching the 
surface the top of the particle is at about 700 K, the surface temperature of the binder. This so happens to 
be very much the temperature for first ignition of HMX (see Fig. 23). There will be however a transition 
delay of the HMX particle to full combustion for which its surface temperature is about 900 K. This will 
be dealt with in more details further on. 


In the case of AP, the ignition temperature is around 650 K and the surface temperature for full 
combustion is around 900 to 1000 K. Immersed in the inert binder with a surface temperature of ~ 1300 K, 
the AP should reach combustion as soon as it is uncovered with no transition delay. 


2.0 Propellant Burning Rate Resulting from Component Rates 


Various approaches have been presented in the 
literature to build the propellant burning rate from the 
components’ own burning rates. One of the first Oxidizer 
comprehensive models for composite propellant 
combustion modeling, Ref. [37], had a picture of AP 
and binder burning in parallel with a partitioning of 
the surface between the two ingredients and a 
surface averaged propellant burning rate. This view is 
also found in Refs. [38,39], with Ref. [40] being a 
complete review for AP-inert binder propellants, 
including references preceding the work of [37]. 
A view of a sequential burning mechanism, in which a 
given path goes through oxidizer particles separated 
by layers of binder, Fig. 39, with as a consequence a 
time averaged propellant burning rate, was first 4 
presented in Ref. [42]. 


Direction of burning 


Figure 39: Time Averaged Propellant Combustion. 


This viewpoint is subscribed to in the present work, as was done in Ref. [43], although the details have 
since then evolved into the option presented herein (which seems to be coherent with the latest view of 
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M. Beckstead [45]. Averaging of the components' burning rates, HMX and active binder, in Ref. [41] 
into the propellant burning rate is obtained by assuming that the components melt and mix at the surface 
and then form an average premixed gas flame. Due to the very small thickness of the melt (to the extent 
that melting occurs) layers (of the order of um) and very short residence times in these layers (tenths of 
ms), it is believed that such premixing of the components should not take place. This is also the conclusion 
of [45]. 


Finally, in the recent work of Ref. [44], a surface average is operated for AP-binder interactions whereas a 
time average HMX-binder approach is adopted within the same mixed oxidizers propellant. 


In a randomly packed arrangement of oxidizer spheres of diameter Dox, with the average height through 
the sphere from a given direction (perpendicular to the surface) hox, and for a volume of 1 cm? on the 
surface by 1 cm in depth, the number of particles intercepted along 1 cm of length being N, one has: 


Nh,, 1 cm?/ Nh, 1 cm? = Eox / (1 - Eox ) (1) 


with h, the average binder height between particles and ox the volume fraction loading in oxidizer. Then it 
comes: 


hy = hox ( 1 - ox ) / Eox (2) 
and due to 
N (hox + hy) = 1 cm, N = ox / hox (3) 


For a propellant burning rate v,, the time to run through 1 cm of depth is 


t 1 / Vb.p Tox H tp N hox/ Vb,ox + N Atox ae N h,/ Vb,b (4) 
At., being the (possible) transition delay to full combustion after the top of the oxidizer particle has 
reached the surface. Then the propellant burning rate, as expressed with the component burning rates, is, 


with Eq. (3) taken into account, 


1 /Vop = Cox! Vb,ox + Cot Atox / hox + ( 1 = Eor) / Vb,b (5) 


Consideration of a sphere being traversed randomly along a given direction leads after some computation 
to hx = Dox (7/47. 


In the case of AP- inert binder propellants with no transition delay, the propellant burning rate is 

1 / Vo,p = Sox / Vo,ap + (1 - Gx) / Vo (6) 
The relation ship between volume fraction loading € and mass fraction loading a being 

i= Si Pi / Pp 7) 
DSi =1,1/ pp = Gox/ par + (1 - Ox) / Pp (8) 
that is for example for an 88 % AP - 12 % PB binder p, = 1.72 g/cm’. 


The mass burning rate of the propellant is 


Mp = Pp Vb,p 
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and Eq. (6) yields 


1 / mp = Qox / Map + (1 - Oox ) / Mp 


lll 


ORGANIZATION 


(9) 


In the case of a propellant loaded with aluminum, it is known that the aluminum particles are ejected from 
the surface [46,47] (Ref. [47] being an extensive review of the processes of aluminum combustion) 
and burn at several hundreds of um from the surface. The view of the combustion of aluminum is 
summarized here, see Figs. 40,41,42. 
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Figure 40: Visualization of Aluminum Behavior (From Ref. [47]). 
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Figure 41: Aluminum Combustion. 
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Figure 42: Behavior of Aluminum (18 %) in a Composite Propellant. 


e Aluminum melts at 930 K (inert binder reaches about 1300 K) (active binder reaches about 700 K) 
and individual particles aggregate in part into clusters. 


¢ Sketches of selected frames of high speed motion pictures, Fig. 40, illustrate a protracted (because the 
protective AlO; coating melts at 2300 K) ignition-agglomeration event in which the hottest 
portion of the accumulate inflames and precipitates the complete inflammation-agglomeration of the 
(already hot) accumulate. 


e The Al,O; protective coating collapses into cap, exposing the Al globule (typically ~ 100 um). 


e Aluminum vaporizes at ~ 3300 K and reacts in nearly-spherical flame with CO, and H,O from 
AP-binder flame. 


¢ The caps give Al,O; particles of a few tens um, the spherical flames give Al,O3; smoke of ~ 1m. 

e Not all aluminum particles go into accumulates (packs of original particles sticking together) 
from pockets within large AP particles. About 2/3 of the original particles escape individually from 
the surface, giving as combustion product essentially 1 um smoke [48]. 

The volume fractions being 

box > Ep > Eal 

Eq. (8) becomes 


1 / Pp = Qox / Pox + Oy / Pp + Oa / Pat (10) 
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that is, for example, for a 70 % AP, 20 % aluminum (pa = 2.7 g/cm°), 10 % PB binder (pẹ = 0.92/cm’), 
Pp = 1.84 g/cm’. 

With respect to the burning rate of the propellant loaded with aluminum it is obtained: 

1 / Vb,p am box / Vb,ox le Nai ha / Vb,al F Nai Atal + Eb / Vb,b (l 1) 


The “burning” rate of aluminum can be considered as infinite since it is ejected from the surface, whereas 
its “transition delay” is the time for the binder to regress through the particle height: 


Ata = ha / Vo,» 
Also, for each component N h = & its volume fraction loading. The propellant burning rate becomes 
1 / Vop =; Cas / Vb,ox + Eb / Vb,b + Eal Von = Cox / Vb,ox T ( 1- Eox) l Vb,b (12) 


that is the relation which would be obtained for a corresponding propellant with no aluminum, with the 
binder filling in for it. 


From an energetics point of view the modeling of the combustion of aluminized propellants should 


include at the surface the heat of fusion of aluminum ~ 95 cal/g of Al. One can see Ref. [49] for such an 
approach. 


3.0 HMX (or RDX) (or HNIW) — Active Binder Propellants 


Both HMX and the active binder have independent burning rates. The resulting propellant burning rate is 
given by Eq. (5), where the transition delay has to be evaluated. 


The emerging HMX particle offers to the external heat flux the surface area of the sphere cap which has 
been exposed by the binder regressing at vy» after the time t from first appearance has evolved: 


S = Dox (Vo t ) 

With a being the part of the sphere which is heated by the superficial flux ® , the temperature rise is 

a ( Pp Cp Jox AT / dt (42/3) (Dox /2 = ® n Dox Vout (13) 
The heat flux received is that of the binder flame when the particle first emerges and it transitions to that 
of the HMX flame when the particle has reached full combustion, a transition formula being assumed to 
be 

D = D+ [ (Ts - Tsp) / (Tsox - Tso ) ] (Pox - Pr ) 

with 

Dox = Pp Vb ( Cg Ts - Cp To - Qs Dox. 

The fraction of the sphere heated by the flux is taken to be, with K finally adjusted at 0.1, 

a = exp( - K Dox / €p* ) 

e,* is the thermal wave thickness for which the temperature is at 90 % of its surface value, that is, 


sufficiently close to it, 
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(T -To)/( Ts - To ) = 0.9 = exp( - e,* Vb» / dp) (14) 


When Dox / e)* — 0 the particle is vanishingly small with respect to the thickness of the layer at about the 
surface temperature, a —1 and the particle is heated in its entirety by the heat flux. If Dox /e)* — œ the 
particle is very large compared to the surface layer, a — 0, it is heated on a vanishingly small part. 
When Dox = e,* , a = 0,9, the particle is immersed in the binder layer at T ~ Ts» and it is almost totally 
heated by the external flux. The transition delay from Eq. (13) is then 


Ato” = ( Pp Cp ox ( Dox / 3 Vo ) [ ( Tsox - Tsp ) ] / ( Box -Do )] In( Pox / ®, ) (15) 


One noteworthy feature of this relation is that the transition delay is proportional to the particle size Dox, 
in such a way that inserted in Eq. (5) it renders the burning rate insensitive to Dox, a fact which is observed 
experimentally for HMX-energetic binder propellants [43]. 


Figures 43 and 44 show two examples of propellant burning rate laws with the corresponding component 
rates. It is seen that the propellant burning rate is intermediate between those of HMX and of the binder 
(its being close to that of the binder at pressures under 100 atm is coincidental). 


o HMX 
O Energetic binder (Heat of explosion = 850 cal/g) 


A 70% HMX [20 um] + 30 % Binder 
—— Computed burning rate 


100 Burning rate (mm/s) 


Pressure (atm) 
E ee oe eS ee 


} CATES EE ees E Ss E 


10 100 1000 


Figure 43: Experimental and Computed Burning Rate of a Nitramine Based Propellant. 
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HMX 
O Energetic binder (Heat of explosion = 850 cal/g) 


A 75 % HMX [40 um (50 %) + 4.7 um (50 %)] 
+ 25 % Binder 
—— Computed burning rate 


Pressure (atm) 
10 100 1000 


Figure 44: Experimental and Computed Burning Rate of a Nitramine Based Propellant. 


The model presented above, and, what is important, the mechanisms it takes into account (that is the 
absence of diffusional interaction between the components and the importance of a transition delay for the 


HMxX particles), is quite representative of the experimental results. Some of the details of the making of 
the burning rate are given in Table 15. 


Table 15: Elements for the Evaluation of the Burning Rate (70 % HMX — 30 % Active Binder) 


% of burning 
time due to 
delay 
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It is seen that at low pressures the transition delay has a strong impact and the burning rate of the 
propellant happens to fall close to that of the binder. At high pressure the propellant burning rate is nearly 
the time average of the component rates (without delay the average at 250 atm is 24 mm/s). 


It is observed that, for a given HMX-energetic binder composition, the burning rate law is locked and that 
there is no way to tailor it (as can be done in AP-inert binder compositions by acting upon the particle 
size). Furthermore the pressure exponent is too high to be acceptable for the motor operation. Attempts to 
act upon the HMX burning rate by the use of additives have not been successful. One possibility of action 
however is with additives specific to double-base propellants, lead and copper salts and carbon black, 
incorporated in the active binder. It has been found on Fig. 12 that a moderate super rate effect can be 
obtained. It is seen on Fig. 45 that such a modified binder associated with HMX, and after optimization of 
the amount and of the size of the additives, yields a burning rate law with moderate pressure exponent and 
initial temperature sensitivity. 


Reference Propellant 
propellant + additives 


V,(70 atm) 5 mm/s 12 mm/s 


Å Burning rate (mm/s) 


100 


| Propellant 
+ additives 
K 
10 

[ Reference 

+ propellant 

|: (65% HMX) 

P (atm) 

1 Sea L 1 a E E | > 
10 70 100 1000 


Figure 45: HMX - Active Binder Propellant with Additives. 


Propellants based on the new oxidizer HNIW and binders made from GAP (glycidyl azide polymer) 
have been considered. GAP has been studied in Refs. [50,51] with the techniques employed for the double 
base propellants, in particular its surface temperature under combustion is similar, around 800 K. In the 
actual binder GAP is associated to a liquid nitrate ester, see Ref. [35], the behavior of which is very close 
to a cool active binder. The modelization has been applied to these propellants (in particular assuming no 
interaction, other than the transition delay for the full burning of the oxidizer), it can be seen to give a 
proper account of the resulting burning rate, Fig. 46. As opposed to HMX or RDX which are balanced to 
give CO, H,O and N>, HNIW is over oxygenated to give 


CHN 2012 > 6 CO +3 HO +6N2+3/2 0, 
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Figure 46: HNIW (60 %) — Binder GAP + Nitrate Ester. 


A possible recombustion with the gases from the flame of the active binder could be conceivable, 
this implying an influence of the oxidizer particle size and a possibility of regulating the burning rate as in 
AP composite propellants. The modelization does not confirm this possibility. Comparing propellants with 
different particle sizes would clarify this point. 


4.0 AP-Inert Binder Propellants 


The burning rate of a propellant based on AP and an inert binder such as HTPB is thought to be described 
by the averaging rule of Eq. (9) from the burning rates of the components. In this case the binder has of 
course no autonomous burning rate. Its regression rate is due to the heat flux from the diffusion flame, as 
depicted on Fig. 38, between O, from the AP flame and the hydrocarbons from the pyrolysis of the binder. 
A description of such a flame is given here. 


A column of O; containing gases is ejected from a particle with an efficient diameter D,,* proportional to 
Dox at a velocity v,, such that m, = Pg Vg. This column is consumed by a lateral diffusion characterized by 
a diffusion velocity 


Vait = AqdD/(D/2) (16) 


D being the local diameter of the column, Ay some constant of order 1 and Ð the diffusion coefficient in 
cm/s. The variation of the column diameter is then given by 


dD=-2 Vaite dt $ dx = v, dt 


D dD =- 4 Aa Ð dx / vg 
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Resulting in 


Dox = 8 Aa Ð Xta/ Vg Xa = Dox m,/8 Aap D (17) 


In this relation Dox is the average diameter of a sphere cut randomly by a plane 
Da =2D,, / v6 


When the lateral diffusion of O, and hydrocarbon gases into each others is purely laminar the diffusion 
coefficient is expressed as: 


D=D)T*/p,pD=DoT* 'M/R 

(with the equation for perfect gases p / p = R T / M being used) and the flame stand-off distance 

X= Dox Mp(R/M)/8 Ag Do T! (18) 
does not depend explicitly upon the pressure. In this limiting case and due to the fact that 

m,~1/x¢, 

it is obtained 

mp ~ 1/ Dox 

independent of the pressure level and directly dependent on the particle size. 

As will be seen later, it has been found that the diffusion flame process might become turbulent at high 
pressures when large differences exist between the mass flow rates emitted from the AP and from the 
binder. This conclusion is also mentioned in reference [40]. A general expression for a turbulent transport 
coefficient is 

pD~pu'l 

with u' the magnitude of the fluctuating gas velocity and its scale. It is then assumed 
£~Dox/2, pul ~ (Mo - Mp ) 


that is the turbulent enhancement is related to the difference between the mass flow rates within the O2 
containing column and in the surrounding gases. Then Eq. (17) becomes 


xa = Dox Mp/ 8 Aa (p D der (19) 
(P B Jat = Ðo T? 'M/R+K (Dox /2) (mo - mp ) (20) 
where K should be of the order of 0.1. 


An extra flame thickness related to the chemical time for the completion of the O,-hydrocarbons reaction 
should be taken into account. It is obtained 


Xi (~V, To, ~MTy,/p ~ m/@)= m,/ Ar p Ags exp ( - Egt/R Te) (21) 


the characteristics being related to the final O- - hydrocarbons flame. 
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The stages of the combustion of an AP propellant are shown on Fig. 47, which is to be seen side by side 
with Fig. 38 for the corresponding values. The temperature profiles are approximated as 


(Tr-T)/(Te- Ts) = exp(- 3 x/ xe) (22) 
with 
qs = 3 des ( Tr = Ts y Xf (23) 
Q, 
T, 
q: 
X 
T Arap 
—— 
nen X; aP Qrar | Asb 
Tsar Tso 
aig AP Binder Q., 
a pe ee Y— = T, 


Figure 47: Stages of an AP — Inert Binder Propellant. 


The mass regression rate of the binder is then given by 

Mp Qep = qs,b = 3 des ( Tr 7 Tsp ) / Xf (24) 
where the constituting elements of this relation have been seen in the chapter on inert binders. It should be 
noticed that, to simplify the description, a uniform mass flow rate m, is taken in the gas phase above the 
binder and the AP flame. The flame height xris then the sum of Eqs. (19 and 21). 

The mass burning rate of the AP, see Eq. (13) of the corresponding chapter, 

Map = [ 3 <@ap> Ags ( Tear - Tsar ) / Qeap lie (25) 
3 <Wap> = p Agar Expl - Egar / R Tear ] 


Now, due to the fact that the AP flame receives a heat flux from the final flame, the flame temperature is 
no longer the adiabatic temperature ( Trap “= 1205 K ). It is given by 


Map Cg ( Trap - Ts ap ) + Map Qe,ap = Map Qe ap + rap (26) 
with, from Eq. (22), applied between X¢,ap and X¢ the heat flux from the main flame into the AP flame 


grap = 3 Ag ( Tr- Trap )/ (Xt - Xap ) (27) 
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Results for the burning rate of AP-HTPB propellants are indicated on Fig. 48, where is seen a strong 
influence of the AP particle size. The model described above, and what is important the mechanisms it 
contains, gives a satisfactory account of these results. 


—— Turbulent diffusion 
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ref 221 (5 microns) 
ref 257 (5 microns) 
ref 262 (5 microns) 
Laminar diffusion 


;ox#*o4 


Burning rate (mm/s) ` 


100; 
10 | 
VA Pressure (atm) 
q Ena a rs 
1 10 100 1000 


80 % AP - 20 % HTPB 


Figure 48: Computed Burning Rates Compared to Experiments. 


For the very small AP particle sizes the final flame is mostly dominated by the chemical process, 
very sensitive to pressure, as is the AP flame. This results in a variation of the burning with a pressure 
exponent close to 1. 


For large particle sizes the final flame is dominated by the diffusion process, which is insensitive to 
pressure. This, combined with the pressure dependent flame of the AP, see Fig. 38, gives rise to a 
propellant burning rate which is moderately sensitive to pressure. It is found in the modeling that 
the contribution of the turbulent diffusion, see Eq. (20), becomes important above about 100 atm. 
However above 300 to 400 atm the burning rate of the propellant undergoes a sharp exponent break that 
the model cannot follow. In some references this exponent break has been assumed to be due to the 
burning rate of the AP which would also increase sharply. It has been seen on Fig. 20 that, when the 
samples are carefully inhibited, this is not the case. 


The strong influence of the AP particle size upon the propellant burning rate, as expressed by the 
modeling, is also shown on Fig. 49. It is seen that at about 1 um there is no further gain in the burning 
rate. This is due to the fact that there will always exist a non vanishing flame stand-off distance related to 
the chemical time for the O2-hydrocarbons reaction. 
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Figure 49: Computed Burning Rate vs AP Particle Size (80 % AP — HTPB). 


On Fig. 50 the burning rate computed from the mechanisms modeled as described above is compared to 
results from the literature [38, 39] for mixed AP particle sizes. The agreement is adequate for the larger 
particle sizes and approximate for the small ones (Now, how representative are the particle sizes 
indicated?). This reveals that the physico-chemical features incorporated in the model of the AP-inert 
binder propellant are probably sound. However, once a model has been “tuned” to represent a set of 
experimental results, as on Fig. 48, it cannot be expected that it will “predict” very accurately other results 
for different values of the parameters. 
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Figure 50: Multi-Modal AP Propellants. 


4-54 RTO-EN-023 


Combustion of Solid Propellants 


It has been seen that by acting upon the AP particle size it is possible to tailor the burning rate of the 
propellant, but that there exist a limiting size below which the effect will be non existent, Fig. 49. It is 
possible to gain further by incorporating metallic additives, Fig. 51, such as of the ferrocenic type 
(which during processing will dissolve into the not yet cross-linked binder for a proper mixing). 
There also exist a polybutadiene binder which incorporates within its chemical structure the ferrocenic 
pattern (to avoid migration problems). Various results, and in particular the similarity of action of a silicon 
binder which produces on the propellant surface a fine structure of SiO, residue, to which it can hardly be 
attributed a catalytic efficiency, indicate that the mode of action of the additive is probably physical 
(rather than catalytic, in the sense of enhancing some chemical reactions). The layer of residue deposited 
on the surface has probably a flame-holding effect, the gases flowing in tortuous paths through this residue 
will react closer to the surface, with the final temperature not being modified, in such a way that an 
enhanced heat flux will act on the surface. This is taken into account in the model by decreasing the 
diffusion flame height, see Eq. (18), by a given factor throughout the pressure range. 
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Figure 51: Additive Modified Propellant. 


Finally, the procedure to describe the combustion of AP composite propellants is applied to a practical 
propellant containing two AP particle sizes and aluminum. It is observed on Fig. 52 that the presence of 
aluminum, which has been seen to burn far away from the surface in such a manner that no significant 
conductive heat flux from its flame zone goes back to the propellant surface, should be taken into account 
by adding a radiative heat flux as measured in Ref. [52] for propellants very close to the ones considered 
here. It is found then that the experimental burning rate law is followed quite well. 
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Figure 52: Burning Rate of an Industrial Propellant. 
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CONCLUSIONS 


A review has been presented of the viewpoint of the authors on the combustion mechanisms of 
components and of solid propellants. Some noteworthy points are stressed here. 


Double-Base Propellants and Active Binders: 


A fairly good knowledge of the processes in the condensed and gas phases has been acquired. The main 
point is the presence of a two-flame system, involving NO. —> NO — N,, collapsing into one flame above 
200 to 300 atm. Specific additives (lead and copper salts and carbon black) have a true chemical 
interaction, that is the enhancement of the NO-carbon reaction, bringing some of the energy normally 
evolved in the second flame closer to the surface. This knowledge carries over to the active binders, 
which are however less prone to super-rate effects because they produce less structured carbon residue. 


Inert Binders: 


Not so many reliable pyrolysis measurements have been performed on actual binders. It has been shown 
(from comparisons for a number of materials) that the pyrolysis characteristics obtained at low heating rate 
(by thermogravimetric analysis or differential scanning calorimetry) should extrapolate and apply under 
combustion conditions. The mechanism of pyrolysis is a thermal breaking ( not influenced by pressure) 
of the cross-links and of the polymer, not affected by additives. 


Ammonium Perchlorate: 


It is believed that the mechanism of combustion is properly known. A large amount, about 70 %, of the 
AP exothermically decomposes in a thin (~ 1 um) condensed phase superficial layer, the remaining 
30 % sublimes into NH; and HC1O, which react in a flame very close (a few um) to the surface. Due to 
the thickness, and the very short residence time associated, of these zones, additives have no true catalytic 
(i.e. chemical) action. 


HMX, RDX and HNIW: 


The combustion of HMX (or RDX) is qualitatively comparable to that of a DB propellant, with the 
occurrence of two flames, involving NO, and N,O/NO. However above 20 atm these flames collapse into 
one and therefore it is not possible to induce super-rate effects with specific additives, as was the case for 
DB propellants. Also, additives which were hoped to accelerate the melting-decomposition of HMX fail to 
act under combustion characteristics. The burning of HNIW, the new caged nitramine, seems to be of the 
same nature, with a condensed phase decomposition due to thermal degradation and probably no 
evaporation phenomenon; also the gas flame is probably very quickly collapsed into a single overall 
flame, with therefore no possibility of additive effect of the double base type. 


AP or HMX Inert or Active Binder (Aluminum) Propellants: 


The approach believed to describe adequately the combustion of composite propellants is a sequential one. 
When following a path through the propellant, it meets successively the combustion of the oxidizer 
particles and of the binder layers. 


In the case of AP-inert binder propellants, the propellant burning rate is an average of the components’ 
rates. However an interaction flame between O, from AP and hydrocarbons from the binder enhances the 
rate of AP and allows the pyrolysis of the binder. The burning rate of the propellant can be quite widely 
tailored by reducing the AP particle size and further by adding metallic compounds which leave an oxide 
residue layer, the action of which is likely to be a flame holding effect, i.e. a physical action. Aluminum 
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incorporated in such propellants of course increases the final temperature by as much as 1000 K, but its 
combustion is so far from the surface that it does not influence much the regression rate of the propellant. 


In the case of HMX — RDX — or HNIW - active binder propellants, both components have autonomous 
burning rates. The resulting propellant burning rate is the average of the two rates, with a further slowing 
down due to the fact that the HMX, or other nitramine, particles, upon reaching the surface, have to 
undergo a transition to full combustion. This delay is of importance up to roughly 150 atm; above it tends 
to become negligible. No additives have been found to act on HMX and additives of the DB type act only 
moderately on the active binder. The tailorability of HMX, or other nitramine, - active binder propellants 
is therefore more limited than for AP propellants. 


Further trends on new ingredients are related in the first place to the need to reduce the vulnerability to 
various aggressions of missile motors employing solid propellants: “cook off’ due to fires, bullet or 
fragment impact, “sympathetic detonation”. Another emerging concern is that of the pollution caused by 
the large boosters of space launchers, using AP - inert binder - aluminum propellants, with production of 
HCI and Al O; smoke. 


These considerations have led to reconsidering ingredients such as ammonium nitrate, NH4NOs, discarded 
previously as being insufficiently energetic, reconsidered and seemingly not considered any longer, and to 
introducing energetic binders such as GAP, glycidyl azide polymer (containing the N; group, inducing an 
exothermic decomposition), briefly mentioned in the text, or oxidizers containing the N-NO> group while 
being less sensitive than HMX or RDX, such as ADN (NHq-N-(NO3),). 


Also considered are caged nitramines, such as CL20 (also called HNIW), which maintains the energetic 
level with a lesser amount of oxidizer, thus reducing the vulnerability, or which allows a higher energetic 
level. A number of new results on this component has been presented in the text. 


It is hoped that the large amount of information gathered and the understanding acquired about the already 
used propellants ingredients will allow an efficient approach to the mastering of the behavior of the new 
ingredients and the corresponding propellants. 
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